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ABSTRACT AND SUMMARY 

Ionospheric propagation effects are a major factor limiting tracking ac- 
curacy of the VHF systems of STADAN (Satellite Tracking and Data Acquisition 
Network). These systems are comprised of Minitrack which is an interferometer 
scheme for precision measurement of direction cosines and the VHF Goddard 
Range and Range-Rate (GRARR) system which measures range, range-rate, and 
angle. The extent of ionospheric perturbations is large or small depending on 
the state of the ionospheric medium and the mode and circumstance of the track- 
ing operation. Under conditions of maximum total integrated electron content 
(N t = 10 18 e/m 2 ) the perturbation in GRARR can amount to about 5 km in range 
and 20 meters/sec in range-rate. These numbers exceed instrumental tracking 
precision by more than two orders of magnitude. In Minitrack, the ionospheric 
perturbations can exceed instrumental precision by a factor of 20. The complex 
behavior of the ionosphere makes these perturbations difficult to predict in 



sufficient detail for the needs of precise tracking. A consequence is the intro- 
duction of a relatively large element of uncertainty into system and mission 
analysis. This applies in particular to evaluating the effectiveness of the cur- 
rently employed ST AD AN ionospheric correction model. The approach of this 
paper to these problems is first to give a parametric analysis of the ionospheric 
perturbations. A special concern in the study is to establish the significance of 
the horizontal gradients within the ionospheric medium. In this regard, journal 
publications are surveyed for observational data pertaining to the level of hori- 
zontal gradients of N t . The interest for Minitrack is that a finite refraction 
prevails at the zenith which, at times, can exceed instrumental errors. A general 
result emerging from the analysis of ionospheric perturbations is the central 
role of N t and its horizontal gradients. The limitation to correction of iono- 
spheric effects on tracking data stems from the wide variability and irregular 
temporal behavior of N t as well as its geographical dependence. An observa- 
tional scheme is sketched for the monitoring of N t at STADAN sites which could 
materially relieve the limitation. 



IONOSPHERIC PERTURBATIONS ON STADAN VHF 
TRACKING ACCURACY 

I. INTRODUCTION 

The ionosphere exerts a degrading effect on the tracking accuracy of the 
STADAN VHF systems. These perturbations, characteristically, are highly 
variable depending on the state of the ionosphere and on the mode and circum- 
stance of the tracking. Accordingly, the ionosphere can be a dominant factor 
in limiting system accuracy or have only a minimal effect. In an overall sys- 
tems view there prevails, due to ionospheric influence, a considerable disparity 
between instrumental design precision and the tracking accuracy. Nonetheless, 
the STADAN systems have a long achievement record for the tracking support 
of scientific satellites. The STADAN hardware performance is well established 
both in basic analytic studies and in field tests. The recent design evaluation 
report on the Goddard Range and Range Rate system (GRARR) gives a detailed 
analysis on instrumental metric accuracies (Ref. 1). In the GEOS-I project 
notable contributions are made toward the improvement and verification 
of STADAN tracking accuracy. Among the several intercomparison system 
studies of this project, the accuracy of Minitrack has been demonstrated to 
approach that of the instrumental design (Ref. 2). The STADAN station positions 
have in addition been determined with increased certainty through the geodetic 
tie-in with the Smithsonian Astrophysical Observatory Standard Earth C-5 model 
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(Ref. 3) . These accomplishments would in themselves appear to confer on the 
VHF GRARR system and on Minitrack the status of a precision tracking system. 
Yet this is denied or at least compromised by ionospheric propagation effects 
(Refs. 4 and 5). 

A correction procedure for ionospheric and tropospheric refraction is used 
in STADAN (Ref. 6). The model is based on the ionospheric data published by 
ESSA in Ionospheric Predictions. Some improvement in tracking data is possible, 
thereby, provided the data corresponds to a period during which the ionosphere 
is, in the mean, adequately represented by the predictions. But this is an over- 
simplified view of ionospheric behavior since the medium undergoes wide and 
irregular variations. One can expect that the corrections of VHF range measure- 
ments in particular will suffer from the gross character of the ionospheric data 
and that a large bias may prevail. An underlying shortcoming is the lack of a 
direct observational method on the ionosphere to validate the data. 

In the precise calculation of Minitrack orbits, data is confined to tracking at 
small zenith angles (less than 20°). This restriction is imposed to minimize re- 
fraction despite a program of ionospheric correction. One finds furthermore in 
reports a general prevailing uncertainty in the analysis of refraction effects 
(Ref. 7). 

Although the ionosphere is a major limiting factor in VHF tracking accuracy 
the problem has not received adequate treatment. There is a notable imbalance 


2 



of effort to the detriment of the ionospheric factor. This deficiency should be 
rectified. There is a need for parametric analysis of the ionospheric effects on 
tracking accuracy and a need dictated by goals of tracking system performance 
for appropriate modeling of the ionospheric medium. We furthermore want to 
understand perturbations as a factor in system problems. To lend clafification 
and provide a background for these matters we endeavor now to sketch the iono- 
spheric-tracking interrelation in several of its varied aspects. 

The bearing of ionospheric perturbations on tracking system performance 
is linked with characteristic detailed behavior of ionospheric structure. Certainly 
the most distressing feature of tracking errors having an ionospheric origin is 
the large magnitude of the errors and their irregular changes. Ionospheric 
patterns are noted for their variability and their unpredictability. A basic param- 
eter for the measure of ionospheric perturbation on tracking is, N , the total 
integrated electron content within a vertical column of unit cross-section extend- 
ing through the ionospheric layer. An example illustrating the diurnal variability 
of N is shown (Figs, la, b, c) in measurements of Garriott, et al (Ref. 8). The 
data is derived from observations on the Faraday rotation of signals of Syncom 
III at the University of Hawaii during the fall of 1964. These show an order of 
magnitude variation in the diurnal excursion of N t with a steep morning rise 
and evening decay. The day-to-day variations in the amplitude of the mid-day 
maximum are about 25% and more. 
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Figure 1(a)-Diurnal variation of electron content at the University of Hawaii for 
September 20-25, 1964, inclusive. [From Garriott, et al., Ref. 8] . 
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Figure 1(b)— Diurnal variation of electron content at the University of Hawaii for 
October 9-14, 1964, inclusive. [From Garriott, et a!.. Ref. 8] . 




LOCAL MEAN TIME 

Figure 1(c)— Diurnal variation of electron content at the University of Hawaii for 
November 1, 3, and 4, 1964. These days were selected to illustrate the large day to day 
and even hour to hour variability sometimes observed. Yet, the nighttime and sunrise 
portions remain nearly identical. Individual data points are also shown. [From Garriott, 
et al.. Ref. 8] . 
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Geographic location has a significant influence on diurnal patterns. Observa- 
tions of Schmerling (Ref. 9) during the summer of 1957 made on bottomside iono- 
spheric profiles show very marked dissimilarity for three different sounding 
locations (Huancayo, Talara, and Panama) in the equatorial zone near a common 
meridian; the diurnal maximum at these stations occurred, respectively, at 
morning, midnight, and evening. 

In addition to the diurnal cycle and geographical location, the ionosphere 
varies widely through the course of the eleven year solar sunspot cycle. The 
ionospheric medium is furthermore subject to irregular and abrupt changes under 
the influence of solar and geomagnetic disturbances. These remarks are to point 
out the complex variability of the ionosphere. Indeed to describe the pattern and 
variability is to list in addition to regular features, hosts of anomalies and 
irregularities (Ref. 10 and 11). 

The result of the varigated ionospheric behavior for tracking systems assumes 
differing forms corresponding to the tracking assignment or the application of the 
tracking data. For one, evaluating the effectiveness of ionospheric corrections 
is very difficult. The designation of a mission tracking error model can also be quite 
perplexing. Where there is a critical tracking interval, as say, in support of a 
spacecraft maneuver, the mission planner is apt to be conservative and take a 
pessimistic perturbation model for the analysis. We emphasize that this tendency - 
of exaggerated error estimate - amounts essentially to a downgrading of the 
tracking system expectation where the demand on performance is especially needed. 
The uncertainties associated with the ionosphere also exert an influence on 
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Minitrack, where there is the aforementioned practice of selective use of 
tracking data restricted to within 20° of the zenith. Though refraction 
errors are thereby minimized it is at the expense of a substantial rejection of 
tracking data. That such restrictions can be excessive appears in the GEOS-I 
tests on Minitrack. This will be described below. 

The refraction correction applied to Minitrack data are based on a horizontally 
stratified ionospheric model. With emphasis on small zenith angles, the refraction 
is generally a small perturbation. The character of the ionospheric data is then of 
less significance. However, one may question the validity of the refraction model 
itself. It is evident from the dependence of ionospheric structure on solar zenith 
angle that the medium is not horizontally stratified. This implies a finite refrac- 
tion for tracking at the zenith. The effect has not as yet been studied for Minitrack. 

Now let us look at how these considerations relate to matters on the analysis 
of Minitrack data. Minitrack is routinely relied upon for tracking and orbit de- 
termination as an integral part of satellite experiments . Physical problems like 
these pertaining to the harmonics of the geopotential field and satellite drag also 
enter in the satellite equations of motion. How the ionosphere can interfere in 
orbital computations is shown in a study of Gooding (Ref. 12) on Ariel 2 (1964). 
Using Minitrack data, residuals between observation and calculated orbital 
position were 1' to 2' arc (~ 0.3 mr to 0.6 mr). It was difficult to explain the 
magnitude of the residuals when the instrumental precision of Minitrack (Ref. 13) 
is 10“ 4 in direction cosine or 20" arc (0.1 mr) at the zenith. 
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The discrepancies were attributed to errors, real or apparent, having 
several possible sources, mainly ionospheric refraction, atmospheric 
drag, and uncertainty in the magnitude of the tesseral harmonics entering in 
orbital calculation. In the effort to evaluate these factors, Gooding, left the 
question of ionospheric refraction substantially unresolved. Significantly this 
was a period hear solar minimum when ionospheric effects are smallest. This 
shows the uncertainty prevailing on the part of an investigator in approaching the 
ionospheric problem. 

As a part of the GEOS-I program, a study was made to validate the accuracy 
of Minitrack measurements (Ref. 2). The test data were tracking observations 
on the GEOS satellite over a five day period, Dec. 31, 1965 to Jan. 5, 1966. This 
was in an early ascending phase of the solar cycle (Wolf sunspot number ~40) 
when the ionospheric electron content is still relatively low. The satellite was 
in a near drag-free orbit (perigee ~1200 km and apogee, ~2500 km). A standard 
orbit based on highly precise optical observations was available for intercom- 
parison with Minitrack test data. The treatment of this data differed in several 
respects from the usual routine of Minitrack. Although accurate orbital calcu- 
lations are limited to data observed within 20° of the zenith, the data acquisition 
in the test was extended to zenith angles of 75°. In addition, ionospheric refraction 
corrections were not applied to the data. The result in general showed a 
gratifying agreement in the residuals between Minitrack data and the optically 
determined standard orbit to within ~ 0.2 x 10 -3 which compares with instru- 
mental precision of 0.1 x 10~ 3 in the direction cosine. 
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This is a substantial corroboration of Minitrack precision. To achieve further 


refinement in Minitrack orbits, the ionospheric refraction must be corrected. 

But to validate Minitrack data we have yet to resolve the significance of hori- 
zontal gradients for refraction. The emphasis on observations near the zenith 
places importance on a general analysis of refraction effects at small zenith angles. 

A principal objective of this report is therefore to investigate the contribution 
of ionospheric horizontal gradients to refraction and also their significance for 
GRARR ranging and range-rate perturbations. As part of this study we review 
ionospheric observational data relevant to these gradients. We describe the 
level of the various ionospheric perturbations on STADAN tracking accuracy. 

The theory of refraction and of Doppler frequency shift in a heterogeneous 
ionospheric medium given herein follows substantially that of Al'pert (Ref. 14) and 
Mityokova and Mityakova (Ref. 15), Tropospheric effects in STADAN tracking 
have been studied by Freeman (Ref. 16) and are not considered here. 

The principal difficulty for making meaningful corrections to ionospheric track- 
ing errors stems from the variability of the medium. In order to cope with this 
problem there is no choice but to monitor suitable parameters of the ionosphere. 

A general result emerging from the analysis on tracking perturbations is the 
central role of N t and its horizontal gradients. This then guides us to what 
should be monitored. We suggest in the concluding section with brief remarks 
how a synchronous orbiting satellite might serve in support of this monitor. 
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We treat in the next part the ionospheric perturbations on range, refraction, 
and range-rate (Doppler frequency shift) in that order. In the succeeding part we 
examine ionospheric data from the literature on horizontal gradients or N t . 

H. IONOSPHERIC TRACKING PERTURBATIONS AT VHF 
Range 

The VHF system of GRARR performs ranging by measurement of the pro- 
pagation delay of sidetones impressed on the carrier signal. In the two two-way 
propagation, the uplink carrier frequency fj is approximately 148 MHz and the 
downlink frequency f 2 approximately 136 MHz. The ionospheric range pertur- 
bation AR c is given by the integral 



( 1 ) 


taken along the ray path S : where c is the speed of light, ds the element of path 
length, s c the path length between station position o and satellite position c, and 
v gl and v g2 the group velocities of the carrier signals. 


Group velocity is determined by the relation (Ref. 22) 


c , Bft 

— = M f 

v g 


( 2 ) 


where ft is the refractive index of the ionosphere. Neglecting the earth's magnetic 
field, ft is 

S (3a, 
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N e is electron density (e/m 3 ) and f is in mHz. To a sufficient approximation at 
VHF 

40.3 N 

/x=l — — - = 1 - 77 N e (3b) 

The numerical values of the parameter 77 are 

uplink Vl = 18.4 x 10” 16 , = 148 mHz 4 

downlink rj 2 = 21.8 x 10 ~ 16 , f 2 = 136 mHz 

From the above equations, we have for the range perturbation 


AR C = 77 


f 

Jo 


N ds 

e 


(5a) 


where (in this section) 77 = 1/2 (rj l + tj 2 ). The integral of Eq. (5) is an integrated 


2 

electron content along the ray path, denoted by n t (e/m ). Thus 


AR = 77 n. (5b) 

C ' t 

The bending of the ray path, being second order in the refractivity, can be 
neglected and n t evaluated along the rectified path. 

z c ^ 700 km. Consider the case of satellite altitudes z c above the ionosphere. 


We write Eq. (5) in the form 


AR c =77 



N 


dr 

COS 4 > 


where 


dr = ds cos 0 

r = radial distance from the earth center 

4> = angle between ray path S and the radial line 


( 5 ) 
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r 0 = earth radius 


r 

C 


The predominant contribution to the integral of Eq, (5) comes in an altitude 
interval about N max (maximum of N e ). Since the geometric factor (cos 9 ) _1 is a 
slowly varying function we can make the approximation of removing the factor 
from under the integral sign and assigning an appropriate mean value correspond- 
ing to a point m: 


V 


AR c =• . 

cos 6 


f"- 


dr 


V N t («.) 

COS 4 > 
r m 


( 7 ) 


The point m is at the median for the electron content along the path 
from o to c ; 6 is the polar angle at the earth center of the point m relative to 
o. We designate m as the ionospheric point and its projection on the earth’s 
surface as the sub-ionospheric point corresponding to the ray path. Thus, AR c 
can be expressed in terms of N t (8 ), the total integrated electron content at 
the sub-ionospheric point. 


In the northern mid-latitude region the altitude z ra is variable, but typically 
about 350 km during daylight hours. Due to uncertainty in determining this 
altitude there is an error introduced in the evaluation of Eq. (7). This we now 
shall consider. From the geometry of Fig. 2 


and 


r.sin 4> = r 0 sin 4> 0 = r m sin <f> m 


( 8 ) 


tan = 


rdd 

dr 


dx 

dr 


( 9 ) 
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c 



Figure 2— Plane of incidence. 
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where x ^is the horizontal curvilinear distance. Using Eqs. (7) and (8) we have 
as the error in AR c from an uncertainty Ar m ( =Az J 


-r,N t 


tan 2 0 A r BN. tan cb 

Mn m t Mn 


cos & r 

m m 


+ 1)- 


Bx cos 


Ar 


( 10 ) 


The relative error is 


or i BN 
- tan 2 <}> + . — — — tan 4> Ar 

~ m _ H.T ' m 


N t Bx 


(10a) 


It is reasonable to take A r m = 50 km since this representative of the uncertainty 
in r m . At zero elevation angle, (10a) is a maximum. Then with z m = 350 km, 

Eq. (7) gives tan $ m = 3. The first term in (10a) is 7%; and the second is 1% for 
an above average logarithmic gradient l/N t * BN t /Bx = 5 x io -4 (km) _1 . 

Bounds of N t . To see how the ionosphere controls the magnitude of AR c we 
need to know the bounds on the variation of N t . There are numerous factors 
influencing N t , the principal ones being the eleven year sunspot cycle and the 
diurnal cycle. There has been observed, in the northern mid-latitude region, a 
very useful empirical relation between the mid-day peak, (N t ) max , of N t and the 
mean sunspot number R. Fig. 3 (Ref. 17) shows the composite of data from a 
number of investigators. There is an approximate linear relation between (N t ) max 
and R for the seasons of the year for R > 40; and for R < 40, (N t ) max levels off 
at minimum of about 1 x io 17 e/m 2 . 
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The maximum of (N t ) max is 10 18 e/m 2 . The observed and predicted R vs time 
is given in Ionospheric Predictions. This dependence is shown in Fig. 4 (Ref. 18). 


The diurnal variation of N t is about an order of magnitude. We can then 
give approximate extreme limits on the variation of N t through the eleven 

Ifi o v ✓ iq . 9 

year solar cycle. These are shown in Table I, with 10 e/m N t £ 10 10 e/m . 


Table I 

Approximate Extreme Variations of N t (e/m 2 ) 
(Northern Mid-latitude Region) 



Diurnal Max. 

Diurnal Min. 

Sunspot Max. 

10 18 

10 17 

Sunspot Min. 

10 17 

10 16 



Figure 3— Sunspot dependence of midday electron con- 
tent. [From Yeh, et al. Ref. 17] . 
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SOLAR CYCLE 20 
BEGINNING OCTOBER 1964 



saiawnN lOdSNns oiuioom 


Figure 4— Predicted and observed sunspot numbers. [From Ref. 18] 




The dependence on zenith angle is obtained by combining Eqs. (7) and (8). 


Thus 


AR 


C 


■n N, (».) 



(ii) 


or 


■fi - 0.88 sin 2 ~ 0 


( 12 ) 


for the mediam altitude z m = 350 km. At the zenith the range AR c is, according 
to Table I, between the approximate bounds 20 meters and 2000 meters. At the 
horizon AR c is 2.7 times that at the zenith. 


This perturbation on AR c is to be compared to the GRARR instrumental 
precision which, at the highest sidetone (for VHF), 20 kHz, is nominally 14 meters 
(Ref. 1). 


Refraction 

In this section the theory of refraction in a heterogeneous ionospheric 
medium (Ref. 14) will be applied to Minitrack. 

It is convenient to describe the refraction with respect to the plane of 
incidence, I , defined by the observer, o, satellite position, c , and the earth 
center. The geometry is shown in Fig. 5. In the presence of horizontal 
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Figure 5— Three dimensional refraction geometry. 
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ionospheric gradients, the ray path S from o to c deviates from the plane I, 

The projection of S onto I is denoted by Sj, while S 0 is the line-of-sight. The 
angles d and <£ in the plane I are as previously defined in connection with 
Fig. 2. The horizontal curvilinar coordinate y measures the small deviations 
of S normal to I , such that x , y , and z form a right-handed system. 

At the observer, the refraction of the zenith angle d 4> 0 measured in I is 
given by 

d 4> 0 =4> n -4> 0 (13) 

where 4> n and <£ 0 are zenith angles respectively of Sj. and S Q . The refraction 
angle of S at o normal to I, that is the angle between S and Sj , is denoted 
by da 0 . The relation to the azimuthal refraction angle a 0 is 


da 


o 


da- 


sin <£„ 


(14) 


Minitrack determines the angular position of a satellite from observation on 
the arrival direction of the satellite beacon (downlink, f t = 136 MHz) signal. A 
crossed pair of interferometers having a east-west (equatorial) and north-south 
(polar) orientations measure direction cosines, cos y and cos /3, respectively. 
The connection between Minitrack observations and the refraction theory is 
through the equations (Fig. 6) 
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( 15 ) 


cos yS = s in 4> Q cos a Q 
cos 7 = sin <p 0 sin a Q 
where a 0 is the azimuthal angle. 

The calculation of the refraction is based on the Fermat Principal that the radio 
phase path <5 between emitter and receiver is minimum, that is, 


8 $ = 


2irf 0 


8 



(16) 


where the index of refraction p. = 1 - 7 jN e (x, y, z). (For Minitrack, = t ? 2 of Eq. 
(4).) The calculation is given in Appendix I where it is shown that 


d^ 0 = “ 


tj r 0 r c sin 4> 0 cos 4 > c 


R 


Jt 0 r 2 cos 3 <t> J 0 


BN e<R -R> dR 


3x 


R 


cos 4 > 


(A22) 


da 0 = 77 


r 


H (R c - R) 

R,. 


dR 


(A26) 


The integrals are taken along the line-of -sight path, S 0 , and R is distance from 
o. These expressions have been given by Mityakova and Mityakova (Ref. 15). . 
The restriction is that 

cos 2 4> » rj N e (17a) 
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V 


(17b) 




« COS 2 4> 0 


and 


V 



(17c) 


The first term on the right side of Eq, (A22) is that associated with a 
spherically stratified medium (Kef. 19, 4). For tracking above the ionosphere, 
z c > 700 km, this term, denoted as d^ 1 ^ , can be approximated as 




V r o r c sin 4>o cos 0 C N t (6 m ) 
R r 2 cos 3 6 

c m 


(18) 


where parameters with subscript m are evaluated at the ionospheric point along 
the ray path. For satellite within the ionosphere N is replaced by the sub- 
satellite electron content, with a small correction to the position of point m. 

The relative error introduced in the approximation of Eq. (18) for an un- 
certainty, Ar m , is 


Ar BN. i 

(2 + 3 t an 2 4> ) - — — + t an <£> A r — — . 

V m/ r m dx N 

m t 


(19) 
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Taking Ar m = 50 km, tan 4 m = 1 and l/N t 3N t /Bx = 5 x 10" 4 (km -1 ), this gives 
errors of 4% and 2.5% for the first and second terms, respectively. 

For small zenith angle, 4 0 < 20° the expression for d<?/ 1 ) may be further 
simplified to 


d *$*>=- 


r]N t tan <t> Q (r Q + z c ) 


r 0 Z c 


( 20 ) 


Thus d<£ 0 vanishes at the zenith as tan 4 0 . 

The maximum of d^ 1 ) is at some altitude generally below 700 km depend- 
ing on the vertical profile of N e , while in the asymptotic limit for large z c , 
d^ 1 ) = r]N t tan 4(/*o • This limit is about an order of magnitude less than the 
maximum value at ionospheric altitudes. 

There is given in Table 2 the perturbation on the cosine of the elevation 
angle, d (sin 4 0 ) (1) = cos 4 0 d^ 1 * (based on Eq. (18)) as a function of satel- 
lite altitude for zenith angles, 4>o = 20° and 45° and N t = 10 18 e/m 2 . The tabu- 
lation can be used in conjunction with Table 1 to estimate variations of the 
refraction effect corresponding to extremes of the 11 year solar cycle and 
diurnal cycle. For an earth orbiting vehicle of, say, 1000 km altitude, the 
following are some general conditions under which the perturbation is apt to 
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Table 2 

- d (sin ) x 10 3 for N t = 10 18 e/m 2 


Z c (km/0 o 

20 ° 

45 ° 

350 

~ 0.4 

,~.2 

750 

~ 0.34 

~ 2.0 

1000 

0.28 

1.6 

2000 

0.17 

1.05 

3000 

0.12 

0.83 

r o 

0.08 

0.61 

2 r 0 

0.06 

0.50 

00 

0.04 

0.39 


exceed 10 ~ 4 . At q> = 45° for daytime and nighttime ionosphere of the sunspot 
maximum and also for the daytime of the solar minimum; while at qp 0 = 20° for 
the daytime ionosphere of the solar maximum and also for the midday ionosphere 
of intermediate sunspot number. This then is descriptive of refractive errors 
in Minitrack due to the first term of Eq. (A22). 


The second term in Eq. (A22), which we denote asd^ 2) , and da Q of Eq. (A26) 
give, respectively, the refractive contributions of the horizontal gradients in and 
normal to the plane of incidence. As these are similar in form it will suffice to 
discuss only one. We choose 


I 

d</>< 2 > = T) f 


5 N. <R„ - R) dr 


Bx 


( 21 )' 


c COS 2 0 

We assign to the geometric factor in Eq. (21) a mean value corresponding, 
approximately, tom', the median point of BN/Bx along the integration path. 
Thus 
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3N e 

Bx 


dr 


( 22 ) 


= 


v(K-K>) 


R cos 2 6 


l S 

I •'r. 


or for z c above the ionosphere 


d 4>i 2) = 


77 (R -R m< ) 3N t 

R e COs2 3X 


(23) 


where R m , is the distance from o along S of the point m' corresponding to the 
altitude z m ,. Note that cos <£ m , is determined by Eq. (7). Thus Eq. (23) relates 
d<^ 2 > to the ionospheric parameter, BN t /Bx and z m , . Consideration of iono- 
spheric behavior shows that the gradient of N t varies extensively, while z m , 
has a more restricted variation. The latter parameter would appear to be 
approximated by z m , since at least for simple ionospheric theory like 
the Chapman model, BN e / Bx ~N e . Accordingly, by Eq. (23) the principal 
ionospheric parameter needed for calculating d^ 2 > (and likewise, da Q ) i S the 
gradient of N t when satellite is above the ionospheric layer. 

At lower altitudes the treatment of Eq. (22) is modified. As a rule, d 0^ 2) is 
not large; the maximum for the strongest observed ionospheric gradients 
is within an order of magnitude of Minitrack precision. Furthermore, for z c de- 
creasing to ionospheric altitudes, dc£< 2 > decreases substantially. In this altitude 
regime aij expression of d<£( 2) corresponding to form of Eq. (23) has as its 
gradient factor Bn t (z c )/Bx and the median point corresponding to R m , shifted. 
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Accuracy requirement on these parameters can be relaxed considerably due to 
the diminished magnitude of d<£< 2) . Accordingly if BN t /Bx and R m , of Eq. (23) 
are known there should not be serious difficulty in treating the lower altitude case. 

The main assumption of the ionosphere leading to these remarks is that the 
horizontal gradients of N e are near maximum in altitudes near z m . To get a better 
look at this, one needs observational data on vertical and horizontal profiles of N e . 
However, such data is scarce, at least in journal publications. We do find top- 
side and bottomside data, but at differing periods of observation. The sample we 
shall discuss gives latitudinal gradients of N g (z) for daylight hours. These are 
the conditions under which the horizontal gradients are generally largest. Observations 
of the bottomside profiles by Wright (Ref. 20) were made during IGY. The ob- 
servations were near the 75° W meridian in the latitude interval 15°N to 50°N. 

The data show mean ionospheric behavior for the month of Aug. 1959 at 1700 
LMT. The latitudinal gradient is negligibly small except for the altitude interval 
of 50 km bordering h 

max 

Extensive synoptic observations of the topside have been made by Alouette I. 
Data on latitudinal variations of N e profiles near the west 75° meridian are re- 
ported by Chan (Ref. 22). Two sets of observations under differing geomagnetic 
conditions are shown in Figs. 7a and 7b; respectively, one for June 7, 1962 at local 
time between 1100 and 1430, and the other for Oct. 31, 1963 at local time between 
1600 and 2200. These data show the maximum level of N e at the geomagnetic 
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Figure 7(a)— Latitudinal variation of electron 
density on June 7, 1963 (Pass, 3432; 1757 
to 1836 UT; Kp = 4" to 3°). [From Chan, 
Ref. 21]. 
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Figure 7(b)— Latitudinal variation of electron 
density on October 31, 1963 (Pass 5427; 
2239 to 2326 UT; Kp = 1"). [From Chan, 
Ref. 21], 
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equator with gradients extending into the midlatitude regions. The horizontal 
gradients of N are seen to be a maximum near h . The median altitude of 
the gradient, z , , is about 100 km above h . If we consider a model with the 
combined features of data of Wright and Chen then we would have z , ~ h „ + 

w m max 

100 km. The example is at least consistent with the notion z m , = z m . 


Criteria on gradient of N t . The evaluation of d<^/ 2) , according to Eq. (23) 
depends on BN t /Bx. The magnitude variations of this ionospheric parameter 
will be the subject of the next section. Here we will establish criteria indicating 
whether the gradients are a relevant factor for perturbation on Minitrack 
observations. For making an apt comparison with Minitrack precision, we 
consider d (sin 4> 0 ) (2) (* cos 4> 0 dcf>^ 2 ^ );that is, 


77 cos <p n (R - R_/) BN, 
d (sin 0 O )< 2 > = — — — -tr- 
ees 2 6, x 


(24) 


This term, in contradistinction to dc^^ 1) , does not vanish at the zenith and further- 
more increases asymptotically with range. 


For a rough numerical estimate of Eq. (24) it is acceptable to write 


d(sin<£ 0 )< 2 > = 2.2 x 10" 5 


c o s 4> 0 
cos 2 4> m , 




) ^ 

Bx 


(25) 


where we have used 77 = 17 2 ; the approximation, R m < /R c = z m , /z. c ; and BN t /Bx 
is in gradient units (G.U. = 10 13 e/m 2 km). Representative cases of interest 
are: <p 0 = 0 ° and 45°; z c = 750 km and > 3r 0 . There is in Table 3, for these 
parameter values the minimum gradient of N t for which d <£< 0 2 > equals the Mini- 
track precision, 10 " 4 . 
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Table 3 

BN t /Bx (G.U.) for d(sin <p 0 ) (2) = 10 -4 


^0 /z c 

750 km 

~3r 0 

0° 

9.0 

4.4 

45° 

6.4 

3,2 


Range -Rate 

In the GRARR system the satellite range-rate measurement is performed 
through a coherent two-way Doppler technique. But the system measures effectively 
the Doppler frequency shift of the uplink carrier signal frequency, f t (Ref. 41). 
Hence, we consider the one-way Doppler frequency shift, f d , which in the non- 
relativistic approximation is for a time dependent ionospheric medium 


_ 1 d 
fj c dt 




p.(j, t ) ds 


Kelso (Ref. 22) has demonstrated that 



M(r c , t) v? e 


c 



(26) 


(27) 


where r c is the unit tangent vector to S at c . The first term on right side will 
not in general vanish at the point-of-closest approach in a three dimensionally 

heterogeneous ionosphere (as it would for a spherically symmetric medium) but 
instead vanishes at another point of the trajectory displaced from the former by a 
finite distance. The effect of the time dependent ionosphere is expressed by the 
second term in Eq. (27). 


30 










Utilizing the relation between the Doppler effect and satellite range-rate, the 
effects embodied in Eq. (27) can be represented as the sum of an unperturbed 
vacuum Doppler term, R c /c, and a term of ionospheric perturbation, AR c /c; that is 



Correspondingly we rewrite Eq. (27) 


(28) 


• • • 

f d <^c - !) ( f c- G o)- 7 c 1 f c ^ J 

~ 7 ~ = + + +— ds 

I, c c c c J at 

1 J o 

where u Q is the line-of-sight unit vector between o and c. From the relation 

R c = u Q ♦ r c we have from Eqs. (28) and (29) 


(29) 


AR = - 1) 7- -T + (f 



(30) 


The effect due to local refractivity at c in the first term on the right is by 
Eq. (3b) .A _ . 

O c - !) V r c =-V N e( r c> *> R c ( 31 > 

where we use the approximation r c ~ u 0 . (77 = rj 1 for formulas related to the 
Doppler measurements, hence to all formulas in this section where 77 appears.) 

The second term of Eq. (30) is due to the effect of the finite refraction angle 
at c. We express this in component form (Fig. 8 ) 


*c ' ( f c - G 0> = v l Xc = ( R C ^ 0 ) dc ^c + y c da c 


(32) 
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where 


X c = angle between f c and u 0 

S* A a 

= projection of r c on -r - u Q 
d 4> c = refraction angle at c between S Q and S t 
da c = refraction angle between S and I-plane at c . 

The Eqs. (30) and (32) show how the components of refraction of the ray path 
at c enter in the formula for AR c . 

The formula for the refraction angles are given in the Appendix B. 


d<p c = - 


rj r Q r c cos 4> 0 sin 4> c 


R 


f 


(N e (r) - N e (? c » dr 
r 2 cos 3 0 


- V 


f C R_ 

J 0 R c 3X 


dR 

COS 4 > 


(A28) 


and 


r c R 3N 

da = - 7 ] I — — — dR (A29) 

J K <>y 

Jq e J 

Inserting the terms of Eqs. (31) and (32) into Eq. (30) and using Eqs. (A19), (A28), 
and (A29) along with the geometric relation (Fig. 9) 

z c = K cos 4> c - R c 4> 0 sin ( 33 ) 
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Figure 9— Velocity components in the 
I -plane. 


we obtain 



re 1 x f C R BN e dR , f 

( ° ^°) J R dx COS cf> +Vc j 
•'0 c ^ •'o 


r 3N. dr ] 

J 


(The Roman numerals will serve to identify the various term of Eq. (34) for 
the discussion to follow in later paragraphs.) 
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It is worth noting that the expression of Eq. (34) can be also derived from 
Eq. (28) without introducing refraction angles (Ref. 15). Rewrite the integral 
of Eq. (28) 


_d_ 

dt 



d f A6(x, y, z) dz 

dt J o cos 4 > 


( 35 ) 


The point c and therefore the ray path are functions of time. Hence the parameters 
under the integral sign x, y and implicitly are time dependent. Thus, using Eq. (8), 


_d_ 

dt 



_d_ 

dt 



jU;{x( t), y (t), z} dz 



r o cos2 ( t ) 
( T o + z ) 2 


(36) 


where the time dependence is explicitly exhibited. By formal application of the 
time differential operation indicated in Eq. (36) and using geometric relations, 
the expression in the parethesis of Eq. (34) follows directly. 

Discussion of AR c . Table 4 lists approximate formulas for estimating the 
terms (H-IV) ofEq. (34) at satellite altitudes, 700 km £ z c £ r Q , and large eleva- 
tion angle. The formula for term V requires z c £ 700 km. There is also given an 
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approximate magnitude level/ the terms can attain under the condition of maximum 
ionospheric electron content and maximum horizontal gradient.* The term values 
are to be compared with the instrumental range-rate precision of 6 cm/sec at a 
sampling interval of one sec (Ref. 1). 

Term I is the local medium effect, which is maximum when at the altitude, 
z c = h max . Term II is a maximum at the top of the ionosphere. This vanishes, 
however, at the zenith and at the point-of-closest approach (<p 0 = 0). Term III 
which is dependent on gradients in the I-plane is commonly an order of magni- 
tude less than term n. It vanishes at <Po ~ 0* Term IV is comparable in mag- 
nitude with III, but need not vanish at the zenith or at the point of closest 
approach. 

At lunar range, the satellite angular velocity is so small that terms I to IV are 
negligible, in which case V can be the dominant term. The latter term is then a 
measure of the time rate of change of electron content along the essentially fixed 
line-of-sight. This may be due to short term fluctuation of the ionosphere or due 
to the diurnal cycle. For the steep diurnal changes of N t and low elevation angles, 
the term can amount to about 50 cm/sec. 

The approximate upper magnitude levels of the various terms of Table 4 

\ 

is one to two orders of magnitude greater than instrumental error. 


*The magnitude of horizontal gradients is treated in the next section. 
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Table 4 

Terms of AR c , Eq. (34) 


Term 

Term Approximation 
(m/sec) 

Upper Magnitude Level 
(m/sec) 

I 

II 

III 

IV 

V 

2 N „ z„ 
cos 9 C 

10" 1 N t r 2 sin 2<p 0 *<p 0 
r 2 cos 3 cp 

m 1 n 

„ 3N t R_ 

2 x 10" 2 “ -<p 0 

Bx COS 2 CO 
* m 

BN f R y_ 

2 x 10- 2 , 1 “ ° 

By R c cos cp m 

5.5 x 10" 2 3N t 
cos cp Bt 

* m 

20 

20 

3 

3 

0.5 

Parameter 

N e 

Nt 

Bx 

BN^ 

By 

BN t 

Bt 

y c 

*c 


R c , R 

c m 

r 

m 

AR C 

Unit 

10 12 e/m 3 

10 17 e/m 2 

10 GU 

10 17 e/m 2 hr 

km/sec 

mr/ sec 

10 3 km 

m/sec 


III. HORIZONTAL GRADIENTS OF N t 

The concern in this section is to assay the magnitude range and variability 
of horizontal gradients of N t . The study is primarily a survey and examination 
of data in the literature. We have in mind a comparison of the observed gradient 
levels with that given in Table 3. A level of 10 GU and greater means that the 
refraction will exceed Minitrack precision for satellite altitude above 700 km at 
the zenith. 
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Observational Methods . Historically, N t was first measured from observa- 
tions on the Faraday rotation in lunar radar experiments (Refs. 23, 24). With the 
advent of space exploration this technique has been largely superceded by obser- 
vation on radio beacon signals of earth satellites. The common methods are to 
observe the Faraday rotation or the Doppler shifted frequency of the transmitted 
signals. To resolve an ambiguity in the determination of N t , the preferred pro- 
cedure in these methods is to employ two harmonically related frequencies, or a 

pair of closely spaced frequencies. The signal observations then give the integrated 
electron content along the moving ray path. For satellite altitudes above the ionosphere 
the data yield essentially N t corresponding to the sub-ionospheric points of the ray 
path. The horizontal gradients thereby obtained will be important in our study. 

For a series of satellite passages, data on N t can be acquired at different 
local times and covering a geographic area about the site. In a limited period of 
observation there is little secular change in the ionosphere and one determines 
a distribution function for N t . Some investigators represent their data as a com- 
posite in contours of constant N t in a plot of latitude vs local time. There is a 
tacit assumption that the diurnal pattern of N t in the region of observation is in- 
dependent of longitude. With this limitation, one can then use the plot to recover 
mean latitudinal and longitudinal gradients as a function of local time and latitude. 


The mean latitudinal gradient between latitudes and X. 2 can be determined 
from the relation 


lat. grad. = 


N t <W 2 ) - N t (\) 


(37) 


38 



or, 


lat. grad. (G. U. ) = 


8.5 (N t (\ 2 ) - N t (\ 1 )) 


X-2 - \ 


(38) 


where N t /k is in units of 10 16 e/m 2 deg. 

A longitudinal gradient of N t can be calculated from the dependence of N t on 
local time. The gradient is then given by 


long. grad. = 


N t (t 2 ) - N t (t x ) 


2 7 t r cos k 

m 


t 2 ~ t l 


(39) 


where T is the period of rotation of the earth; or 


long. grad. (G. U. ) = 


5.6 N t (t 2 ) - N t (t,) 


cos k 


t 2 ~ t l 


(40) 


where N t /t is expressed in units of 10 17 e/m 2 hr. We note however that, 
contrary to the assumption stated above, the diurnal pattern is dependent on 
longitude. The ionosphere has consequently an asymmetry which contributes a 
bias to the long. grad. This effect is smoothed out in the composite data. 

Observations on the signals of satellites in a geostationary orbit gives a con- 
tinuous record of the diurnal pattern of N t at a fixed sub -ionospheric point. One 
can also use Eq. (40) to derive the long. grad, contributed by the dependence of 
N t on local time. 
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The deficiency in connection with using Eq. (40) implies that the maximum 
long. grad, will be underestimated. But since the lat. grad, tends to attain higher 
levels than the long. grad, the former is of more interest to us. 

Scale Size of N t Gradients. The records of Faraday rotations and Doppler 
shift observations commonly exhibit irregularities which are attributable to varia- 
tions in the integrated electron content along the changing direction of the ray 
path. The smoothed record gives regular ionospheric gradients that are repre- 
sentative of geographic variations. The irregular deviations from the smoothed 
record are evidence of an ionospheric structure, termed in the literature as 
"large scale irregularities," Both the regular gradients and irregularities are 
to be examined for their significance to VHF refraction. 

The horizontal scale of the irregularities vary from a few kilometers to many 
hundreds of kilometers (Ref. 25). The vertical dimension is typically about one- 
fifth of the horizontal extent. Electron densities are observed to deviate from the 
mean background by anywhere from several percent to a small fraction of a percent. 

There is a minimum size to the irregularities for which the application of refraction 
theory (Fermat's principle) is admissable. The scale size, L, of the irregular region 
must be greater than the first Fresnel zone. If R x and R 2 are distances of the 
irregularity from the satellite and ground observer, respectively, and X is the 
wavelength of the signal then 
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Taking Rj - R 2 = 400 km and \(VHF) = 2m, then L > 1 km, which is about 
the minimum size associated with large scale irregularities. 

Regular Gradients 

Equatorial Region. A few determinations of the latitudinal dependence of 
N were made by Olatunji (Ref. 26) from Faraday rotation records during the 
passages of Transit 4A and 4B for the period May 1962 to January 1963. The 
station location was Abadan (3.9 °E, 7.5 °N, Magnetic latitude 3 °S). The pre-dawn 

observations (Fig. 10a) show the lat. grad, is nil ( < 1 G. U.) , but for the morning 
hours a maximum of N t (Fig. 10b) prevails at 8° to 9°N latitude with steep gradi- 
ents of 5 to 25 G. U. to the north and south of the geomagnetic equator. The 
morning maximum of N t gives way to an afternoon depression in N t (Fig. 10c) 
with lat. grad, high of 25 to 35 G. U. We mention that the approximate asymmetry 
of the minimum near the geomagnetic equator is akin to the "equatorial anomaly," 
wherein the latitudinal dependence of N max (the maximum of the vertical N e pro- 
file) also exhibit a minimum near the geomagnetic equator. At midnight the 
latitudinal variation of N is reverting (Fig. 10c) to a pre-dawn like profile, 
but the lat. grad, remains appreciable, as high as 6 to 20 G. U. 

Tyagi (Ref. 27) measured N from the satellite S-66 for a period during 
the minimum of the sunspot cycle, Oct. 1964 to Sept. 1966. The receiving sta- 
tion at Delhi (28.63 °N, 77.22 °E) recorded data over a wide latitude range. The 
daytime observations (Fig. 11) covering the geographic latitudes 10 °N to 10 °S 
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show lat. grad, of 2 to 20 G. U. The nighttime observations covered the latitude 
interval 16 °N to 40 °N. The lat. grad, were low, < 1 G. U. due to the low level 
of N t at night. (The logarithmic gradients of N t in units of inverse degrees 
are shown in parentheses in Fig. 11.) Tyagi also gives a summary on the 
diurnal variation of N t . These data imply a long. grad. , ~ 3 G. U. 
for daylight hours during the equinox. 

Observations on N in equatorial region during 1964 (sunspot minimum) were 
made by Rufenach et al (Ref. 28). The observations of signals of Transit 4- A 
involved a cooperative effort of ground station at Bangkok, Singapore, and Hong 
Kong. The long. grad, (deduced from diurnal variations) of N t were low, about 
2 G. U. and less. The data on latitudinal variations of N t for May- June 1964 
(Fig. 12a) and Aug. - Sept. 1964 (Fig. 12b) show large lat. grad, during.the late 
morning and early afternoon of 10 G. U. to 30 G. U., but near the magnetic equa- 
tor, where N t is an extremum, the gradients are nil, 

Basu et al (Ref. 29) have measured the latitudinal variations of N t for quiet 
and disturbed geomagnetic conditions. The data is derived from signals of the 
beacon satellite S-66 at Calcutta, 23°N (dip 32°N). The observations on quiet 

days exhibit an interesting regularity (Fig. 13 a). The curves of N t all peak 
around 23° north latitude and decrease toward both the geomagnetic equator 
and higher latitudes. The influence of seasonal effects and solar 
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Figure 1 1— Latitudinal dependence of N r [From Tyagi, Ref. 27] . 
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Magnetic Dip Angle, degrees 

Figure 12(a) “Latitudinal variation of electron content: three-station 
analysis for May and June 1964. [Rufenach, et al., Ref. 28] . 




Figure 12(b)— Latitudinal variation of elec- 
tron content: two-station analysis for August 
and September 1964. [From Rufenach, 
et al.. Ref. 28] , 
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activity appears as a relative displacement of the curves and a systematic change 
in slope. Toward the magnetic equator the observed slopes are 25 to 50 G. U. 
Fig. 13b compares the noontime latitudinal variation of N t for magnetically 
quiet and disturbed days. The significant changes during disturbed days is a 
decrease and flattening of the slope toward the lower latitudes, while at higher 
latitudes the lat. grad, are general high, about 20 to 50 G. U. 



SUtlONOS DIP (,*H) ► 

Figure 13(a)— Variation of electron content 
with latitude under magnetically quiet con- 
ditions in different seasons. [From Basu, 
et al. Ref. 29] . 
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SUilOMOS DIP ( # h) ► 

Figure 13(b)— Comparison of latitude vari- 
ation of electron content under mag- 
netically quiet and active conditions in 
(a) summer (July 1966), (b) fall (October 
1966), and (c) winter (January 1967). 
The solid line curves represent active con- 
ditions, and broken line curves represent 
quiet conditions. [From Basu, et al. Ref. 
29]. 
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The diurnal variations of N t (Fig. la, b,c) observed by Garriott, etalduringthe 
relatively low period of the sunspot cycle, Sept.-Nov. 1964, at Hawaii from signals of 
Syncom IE show a rapid morning rise and evening set (Ref. 1). The day to day 
variations in these rates are themselves highly variable and correspond to long, 
grad, of 1 to 6 G. U. The long. grad, during the night are nil. Coordinates of the 
sub-ionospheric point are 19.7°N, 159. 8°W. 

Mid- latitude Region, de Mendonca (Ref. 30) calculated N t from measure- 
ments of the differential Doppler frequency shift of Transit 2-A. Observations 
were made at Stanford University during the months July-Oct. 1960, a period of 
high solar activity. The determinations of N are shown as contours of constant 
N t in two plots of geographical latitude vs local time, one for July-Aug. 1960 
(Fig. 14a) and the second for Sept. -Oct. 1960 (Fig. 14b). The difference in the 
data of the two figures for their common time interval 1200 to 1600 LT is due to 
the seasonal change in the ionosphere between the periods. The long. grad, 
and lat. grad, are determined as prescribed by Eqs. (38, 40). The gradients 
are seen to vary over a wide range of magnitude and direction. During midday 
the direction is predominantly latitudinal and during the pre-dawn and post-sunset 
hours predominantly longitudinal. 

For the period July-Aug. 1960 the midday lat. grad, varies from about 8 to 
40 G. U. An inspection of the long. grad, at the 30° latitude indicates about 8 G. U. 
at sunset and 6 G. U. at sunrise. During most of the daylight hours the gradients 
of N t are from 5 to 10 G.U., but at night the level drops to 1 to 2 G. U. and less. 
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In the period Sept. -Oct. 1960, the midday level of N t is about 50% higher 
than during the preceding summer months. The midday lat. grad, is uniformly 
high at about 20 G. U. The gradients after sunset also have a relatively high 
long. grad, of 5 to 10 G. U. 

The Joint Satellite Studies Group (Ref. 31) made observations on N t in the 
European sector during October 1963 and during October 1964. The data was 
reduced from Faraday rotations of signals from Transit 4A and Explorer 22. 

The program examined geographical distribution of electron content under differ- 
ing geomagnetic conditions. A unique feature of this effort was simultaneous 
measurement of N from several tracking stations. An example given in Fig. 15a 
shows electron content variations along three sub-ionospheric traces. The Studies 
Group utilized this type data and other ionospheric sounding data to construct tri- 
dimensional (relief) representations showing electron content as a function of 
geographical location. Two such representations for differing geomagnetic con- 
ditions are shown: Fig. 15b, for the quiet magnetic state has gradients of about 
10 G. U. radiating from an area of minimum N ; Fig. 15c, for the disturbed 
magnetic state has gradients of 10 to 20 G. U. having a north-east to south-west 
orientation. These are relatively high gradients for a period of near minimum 
solar activity. 

Determination of the latitude dependence of N t have been made by Weise and 
Yeh(Ref. 32) at the University of Illinois for the period Oct. 1964 to March 1965. 
Observations were made of the Faraday effect on signals of Explorer 22. The 
horizontal gradients of N were consistently very low 1 2 G. U. 
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Figure 14(b)— Smoothed curves of constant / N dh derived from measurements on quiet 
days during September and October 1960. Curve A represents sunset at sea 
level and curve B sunset at height h = 400 km for September 12, 1960. [From 
de Mendonca, Ref. 30] . 


Auroral Region. Liszka (Ref. 34) determined N at Kiruna, Sweden (67.8°N, 
20.3°E) from signals of the satellite S-66. His observation covers the period, 
Nov. 1964 and the Winter 1964-65. Detailed records of N for a series of satel- 
lite passes during Nov. 1964 show lat. grad, of 2 G. U. before local noon and 
after midnight. A representation of observations during the winter are shown in 
a plot of local time vs geographic latitude (Fig. 16). The observed data were 
within ±1 hour of the Kiruna meridian. The diurnal cycle is covered for the 
latitudes 55°N to 72 °N. The notable feature in the contours of constant N t is the 
existence of the minimum which dominates the pattern in the nighttime over an 
extended latitude interval. During the morning hours (~ 0700 to 1100 1. t.) the 
gradients are predominantly longitudinal, about 1 to 3 G. U. in magnitude. In 
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Date id. 10- 64 
E = 17 h 

K = a 

Figure 15(a)-Tracing of the sun-ionospheric lines in 
terms of the geographical latitude and longitude; the 
figures give the total electron content in 10 16 e/m 2 . 
[From Bertin, et al.. Ref. 31] . 
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Figure 15(b)— Tridimensional representation of 
the distribution of the total electron content in 
10 16 el/m 2 . [From Bertin, et al„ Ref. 31] . 
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Figure 15(c)— Tridimentfonal representation of the 
distribution of the total electron content in 10 16 
el/m 2 . [From Bertin, et al.. Ref. 31] . 
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Figure 16— Latitidunal and diurnal variations of 
the ionospheric electron content near the auroral 
zone during winter 1964-65. Units are 10 16 
e/m 2 . [From Liszka, Ref. 34] . 

the late afternoon from 1200 to 1700 1. t. the pattern is irregular with gradients 
of 1 to 5 G. U. The low level of gradients is due to the very low range of N t 
(< 5 x 10 16 ) and characteristic of the minimum in the solar cycle. On the 
otherhand, the logarithmic derivative of N t is large, so that during the more 
active phase of the solar cycle one would look for an appreciably higher gradient 
of N t . 

Rai and Hook (Ref. 33) measured the latitudinal dependence at N t at College, 
Alaska during the winter 1963-64. The signals were from satellite 1963-38c. 

The gradients had a generally low level of about 3 G. U. and less. This again 
corresponded to the low level of N t which was < 1 x 10 17 e/m 2 . 

Large Scale Irregularities of N t , 

Mid- latitude Region. Merrill et al (Ref. 35) analyzed irregularities in the 
Faraday rotation records of the 20 mHz signal of Sputnik. The observations were 
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made at Boulder for the period Sept. 1958 to March 1960, This was an epoch of 
high solar activity for which the integrated electron content N. had a diurnal 
range from a nighttime minimum of 1 x 10 17 e/m 2 to a midday maximum of 
9 x 10 17 e/m 2 (Ref. 36). The horizontal dimensions of the observed irregularities 
were typically about 150 km to 200 km. The excursion within the irregularities 
from the mean background of N t was about 1% to 2.5%. Their gradients in the 
nighttime were about 1 G. U. and the mean for the daytime was 3 G. U., while 
the midday maximum was 6 G. U. 

Titheridge (Ref. 25) analyzed the Faraday fading records of the 20 mHz 
transmissions of Explorer for the period Nov. 1960 to Aug. 1961. The location 
was near Auckland, New Zealand (lat. 30,01°S, long. 174.98°E). Some 700 ir- 
regularities was observed varying in horizontal extension from 5 km to 500 km. 
The data, summarized in Fig. 17, gives the gradients (c/s = 1 G. U.) vs local 
time. The "a" plot is for summer and the "b" plot for winter. In general the 
gradients are very low, less than 1.0 G. U. It should be noted that the level of 
N t was very low. The midday peak is 2.5x10 17 e/m 2 in the summer and 1.5 x 10 17 
e/m 2 in the winter months. 

Roger (Ref. 37) deduced gradients of large scale irregularities from obser- 
vations on the 108 mHz signal of Vanguard at Jodrell Banks. The period was 
Nov. 1959 to Feb. 1960 and July-Aug. 1960. His data show average gradients of 
2.5 G. U. The observations of Merrill, et al, and Roger made at similar times 
in mid-latitude regions agree on the level of the gradients. 
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electron content of the ionosphere. The circles show the variation in the mean value of C/S. (a) Summer; (b) Winter. [From Titheridge, 













Bhonsle (Ref. 38) obtained data on gradients of irregularities for a period of 
low solar activity, Oct. 1964 to Oct. 1965. The deductions are based on the 
Doppler shift of beacon signals of S-66. The stations was at Stanford. The 
average gradient was 1 G. U. 

IV. CONCLUDING REMARKS 

Horizontal Gradients of N t . 

Observational data on the regular gradients of N t have wide variations 
that show, like N t , the influence of time of day, geographical location, phase 
of the sunspot cycle, and geomagnetic disturbances. The largest gradients of 
N t prevail during daylight hours and have predominantly a latitudinal direction 
The general trend of the midday latitudinal profile of N is an irregular decrease 
from the equatorial region to the auroral regions. A common feature of this pattern 
is a depression of N at the geomagnetic equator. 

During the period of lowest sunspot number the latitudinal gradients in the 
equatorial region may attain a level as high as 10 G. U. and during the phase of 
high solar activity 60 G. U. has been observed. For the mid-latitude region the 
observational data indicate latitudinal gradients approaching the levels of the 
equatorial region. At sunrise and sunset, gradients are predominantly longi- 
tudinally directed, and attaining levels of 10 G. U. for the solar maximum. 

In the near auroral region the limited observational data show gradients of 
several G. U. during aperiod of low sunspot number. The region however has strong 
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irregularities and correspondingly large logarithmic gradients of N t , For high 
sunspot number the gradients might attain a more significant level. 

The gradients of N stemming from large scale irregularities are substantially 
less, by about an order of magnitude, than those of regular gradients and show 
the influence of sunspot number. Most observations on the irregularities have 
been confined to the mid-latitude region. During high sunspot number, gradients 
are typically 2.6 G. U. during the day with a maximum of about 6 G. U., while 
at lower sunspot number the magnitudes are generally negligible. 

Refraction 

The effects of a three dimensionally heterogeneous ionosphere on refraction 
formulas is to add to the known terms of a spherically layered medium other 
terms with horizontal gradients. The familiar terms (associated with spherical 
symmetry) attainthe largest magnitude, but they also vanish at points of sym- 
metry, namely at the zenith for refraction and at the satellite point-of-closest 
approach in range-rate tracking. On the tracking interval where these terms 
are accordingly minimal, the horizontal gradients of the ionosphere distinctively 
can contribute a predominant effect, sufficient to warrant consideration in the 
analysis of tracking data. For Minitrack, we also mention that at relatively 
large zenith angles (say, 45°) and large range (several earth radii) the refraction 
terms of horizontal gradients in N t cannot only exceed the limits of instrumental 
precision, but can as well attain a general level of comparability with the term. 
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The latter remark on Minitrack should be qualified because instrumental 
error is a function of signal strength. The radiation level of satellite signals 
is of course power limited. Then for a fixed power level there is a range beyond 
which instrumental errors must increase. Tracking precision can deteriorate 
rapidly with increased satellite range (Ref. 39). This effect eventually over- 
whelms ionospheric perturbation. As radiated power levels of a small fraction 
of a watt suffice for near earth orbit it would appear that an increased power 
level is possible, thus extending the high precision to a range of an earth radius 
and more. It is in this eventuality that refraction due to horizontal gradients is 
particularly relevant to Minitrack. 

Ionospheric Monitoring and Modeling 

We have in the main body of this report given a framework for the analysis 
of ionospheric perturbation on VHF tracking and used observational data to study 
the magnitude and variability of the perturbations. In order though to apply the 
analysis to tracking data, a description of the ionospheric state is needed cor- 
responding to the time and site of the tracking. But the wide and irregular 
variability of the ionosphere dictate that some monitoring of the medium be con- 
ducted in conjunction with the tracking events. We will describe in rudimentary 
form a method for such ionospheric monitoring and modeling. This method, to 
be outlined in the immediate following paragraphs, raises questions concerning 
errors associated with the model or concerning tracking requirements satisfied 
by the method. These, are rather ramified matters. A close statement here is 
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hardly possible. Abundant observational data over extended periods and experi- 
ence with the method are likely to be best guides to answer the questions* 

The method we have in mind is one giving a description of spatial and 
temporal variations of the three dimensionally inhomogeneous ionosphere for 
the individual ST AD AN sites, A local vertical profile N e (h) is to be constructed 
from observed bottomside profile N b (h) and a calculated topside profile N ea (h). 

A model of horizontal gradients in this profile then complete the three dimen- 
sional structures. 

The monitor makes two types of standard ionospheric observations in the 
neighborhood of each tracking site. One set of observations in on the bottomside 
profile. The other set is on N t . The feasibility of the later observation is predicated 
on the expectation that there will in the near future be within the view of each site 
a geostationary satellite whose signals can be utilized for monitoring N t . The 
increased applications role of these satellites makes the expectation appear to be 
a reasonable one. 

With these observations a complete vertical profile is constructed as follows. 

The N e>b (h) gives N max , h max> and the bottomside electron content N b by the formula 

J p“ma x 

N b (h)dh (42) 

o 

For the topside, we assume a model of constant scale parameter, H a , for the 
Chapman profile 
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N =N exp — 'Si 

a , e max .2 


h -K 

2H 


- e 


h - h 

m 

2H 


(43) 


Although the parameter H a is known to be a function of altitude, the profile can 
be adequately described for the purposes of a refraction model by a constant, 
independent of height. 

The topside electron content N a is given by the integral 

N e,a ( h ) dh (44) 



For the Chapman model (Ref. 40) 


N a = 2.13H a N max 


(45) 


Then the relation 


N t = N a + N b (46) 

combined with Eq. (45) gives the topside parameter H a . Thus a vertical profile is 
determined for one location in the site area. 


Note that the bottomside profile and N t which form the bases for determining 
the complete profile, will generally be observed at geographical points having a 
finite separation. In the construction of N e (h) there was the tacit assumption that 
the observations can be identified with a common location. This involves an 
assumption on horizontal gradients, to be considered next. 
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The desired circumstance for monitoring horizontal gradients in the neighbor- 
hood of a given site is to have observations on N e b (h) and N t from pairs of 
stations having an appropriate latitudinal and longitudinal separation. Where this 
is not the case we must look to other means. 

The latitudinal gradients can be relatively large hence their determination 
requires care. A fortunate feature is the location of a number of STADAN sites 
and ionospheric stations near the 75°W meridian. This is favorable for cooperative 
observation and determination of latitudinal gradient. The other STADAN sites can 
be expected to present unique problems according to the disposition (or lack thereof) 
of observing stations in their region. 

Of particular interest for the study of latitudinal gradients are the polar 
orbiting topside sounders. The data already obtained provides a considerable 
basis for developing a model on latitudinal gradients. 

For longitudinal gradients, where stations are not suitable located, there are 
other means as well. As previously stated longitudinal gradients can be attributed 
to the diurnal process and asymmetry in the ionospheric (or geomagnetic) field. 

The part of the longitudinal gradient due to the diurnal effect can be determined 
from the local measurement of the pattern as given by Eq. (40) . To determine 
the contribution of this asymmetry, the global data in Ionospheric Predictions, 
published monthly by ESSA would be a source of data worth studying. 
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The effectiveness with which the general scheme of model and monitoring 
the ionosphere as just described can be used in analysis of ionospheric perturba- 
tion is enhanced by the direct observations on N which is the central parameter. 
A limitation on the model rests with uncertainty in the determination of horizontal 
gradient of N t . 
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APPENDIX A 


Refraction in a Heterogeneous Ionospheric Medium* 


Writing the line element ds in the Fermat principle as 

ds - y&x? + dy 2 + dz 2 (Al) 

or since r = r Q + z and dx = rd0 


ds = 




dr 


one derives from Eq. (16) the familiar Lagrangian equations 


(A2) 





2 B fj. 

By 


0 


(A4) 


*The analysis of this appendix follows closely that of Al'pert (Ref. 14). 
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Multiplying by dr and integrating we have, using Eqs. (9) and (A2) 


/x r sin 4> 


jjL r sin 4> 


iR I 


COS 2 4 > 



(A5) 




dy 

ds 




(A6) 


When 3 yx/3 6 = 0, Eq. (A5) becomes Snell’s law for spherical symmetry, viz., 


/x r sin 0 - (/x r sin cf>) Q - 0 


(AT) 


For 


3/x _ 


0 , 




dy 

ds 



= 0 


(A 8) 


and the ray path is confined to the plane of incidence. 

The refraction will be calculated under the condition of Eqs. (17a, b and c); 
that is. 


d 4 > « 1 
da « 1 


(A9) 

(A10) 
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Eq. (A5) can be written to terms in first order as 


/x r sin <t> - (ji r sin <f>) 0 = 



(All) 


since dy/ds = da « l and where we have the relation for the refr activity 


/* ~ 1 = = - r, N e 


(A12) 


From Eqs. (All) and (9) we have 


dd 

r 

dr 


tan = 



(A12) 


where the refractivity at r = r Q is taken to be zero and the refracted angle 4> at o 
is ^ . 

Using Eq. (13), we have to terms of the first order under the approximations 
(17) and (A9) that 


r dB _ r 0 sin $0 8 r o r2 COS A) t V N e r 0 1-2 sin ^0 

n ( r > + n 3 (r) + 3 (r) 


jliL f ds 

fi 3 (r) Jo 30 


(A 13) 
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where 


H (r) = /r 2 - r 2 sin 2 c£ 0 


(A14) 


Now multiplying by dr/r and integrating from o to c, we get 



d6 = r Q sin <£ 0 



dr 

r n (r) 


+ r o cos d< ^o 


r° rdr 
Jr 0 n 3 (r) 


+ T? r 0 sin 4> 0 



N rdr 

e 

n 3 (r) 


- V 


r c rdr 

Jr 0 « 3 (r) 




(A15) 


As the path of integration can to first approximation be taken along the line-of- 
sight, we obtain from Eqs. (8), (9), and (A14) that 



J ” rc r Q sin <p Q dr 

' . r n (O 


(A16) 


and 


0 (r) = r cos 4> 


(A17) 


Hence the first term on the right side of Eq. (A15) cancels the left side giving 
for the refraction 
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where dR, the rectilinear line element, has replaced ds. 


The denominator of (A18) is evaluated using Eq. (A14) , 


vA18) 


f r ° d r T ° 

2 r 3 , = 2 rdr 2 ^ n3/2 = < r o cos ^o)' 1 - (r. cos )" 1 ( Al9 ) 

J r 0 r co s 3 <t> Jt 0 (r 2 - r 2 sin 2 4> 0 ) 3/2 


The second term in the numeratpr of Eq. (A18) may be simplified by inter- 
changing the order of integration. Thus 


f'-^ — ( f 

Jr 0 r 2 COS 3 0 ^J. 0 


3N 


dR 


r 


BN 

W 


dR 


A 

f 

k J r r z CO S 


OS ■ 4> 


■f 


dR |f- ^(r cos <£)• 1 - (r c cos 0 c )" 1 > 


(A20) 


From the geometry of Fig. 5 we have for the distance R from o along S Q 

R = t cos 4> - r Q cos cf> 0 (A21) 
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Using Eqs. (7), (8), (A20), (A19), and (A21) one can rewrite Eq. (18) as 


d0 o = ~ 


V r Q r c sin 4> Q cos 4> c 



N dr 

e 

— + v 

r 2 cos 3 4> 



(R c - R) ^N e 

R c a* 


■ — (A22) 

COS 4 > 


To calculate da Q we return to Eq. (A6). We mention that Al'pert gives in 
place of Eq. (A6) a form that is symmetrical with Eq. (A5). Our procedure 
simplifies the calculation. 

Since da Q = (dy/ds) 0 to first order terms, we have integrating Eq, (A6) from 
o to c along the line -of -sight 



dy - da Q 



(A23) 


The first term on left side vanishes since the end-points of ray S are on the 
plane I. Hence we have, interchanging the order of integration. 


da Q = 




dR 

~ 


(A24) 
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Then to first order terms in the refractivity 


da 


o 



Finally, 


da Q = 7 } 



3N e (R - R) 

— £ — £ dR 

R c 


(A25) 


(A26) 


The formulas in Eqs. (A22) and (A25) for d<f> Q and da Q are due Milyokova and 
Mityakova (Ref. 15). 

The refraction of the azimuthal angle d a 0 is 


da. 


da Q = 


sin 4> 0 


(A27) 


In connection with the Doppler effect it is useful to have the refraction 
angles at the satellite. We have similar to the derivations leading to Eqs. 
(A21) and (A25) that 
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d 4> c j? + 


1 r 0 r c sin <£ e cos <f> 0 
R c 



(N e (r) - N e (r c ) dr) 
r 2 cos 3 4> 


dr 


~ E ° R dR 

0 R c cos 4> 


(A28) 


and 


d cl c = - 77 



R 

R c 


BN e 

3y 


dR 


(A29) 
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